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^ '. ABSTRACT 

The X-ray spectra of Gamma-Ray Bursts can generally be described by an absorbed 
I— I, power law. The landmark discovery of thermal X-ray emission in addition to the power 

■ law in the unusual GRB 060218, followed by a similar discovery in GRB 100316D, 
HH ■ showed that during the first thousand seconds after trigger the soft X-ray spectra can 

. , be complex. Both the origin and prevalence of such spectral components still evade un- 

i-S^ ■ derstanding, particularly after the discovery of thermal X-ray emission in the classical 

GRB 090618. Possibly most importantly, these three objects are aU associated with op- 
Q . tical supernovae, begging the question of whether the thermal X-ray components could 

be a result of the GRB-SN connection, possibly in the shock breakout. We therefore 
\ performed a search for blackbody components in the early Swift X-ray spectra of 11 

■ GRBs that have or may have associated optical supernovae, accurately recovering the 

thermal components reported in the literature for GRBs 060218, 090618 and 100316D. 

, We present the discovery of a cooling blackbody in GRB 101219B/SN2010ma, and in 

^ ■ four further GRB-SNe we find an improvement in the fit with a blackbody which we 

\ deem possible blackbody candidates due to case-specific caveats. All the possible new 

, blackbody components we report lie at the high end of the luminosity and radius 

Tij" ■ distribution. GRB 101219B appears to bridge the gap between the low-luminosity and 

^-H I the classical GRB-SNe with thermal emission, and following the blackbody evolution 

. we derive an expansion velocity for this source of order 0.4c. We discuss potential 

■ origins for the thermal X-ray emission in our sample, including a cocoon model which 
^SJ , we find can accommodate the more extreme physical parameters implied by many of 
T-H ■ our model fits. 
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1 INTRODUCTION 2010; Nakar & Sari 2012), first suggested by Colgate (1974). 

T, /^T^T^ GRB 060218 displayed indisputable spectroscopic evidence 

The maionty or Gamma-Ray Bursts (GRBs), those belong- „ , ,■ 1 1 ,1 ^t-.- , ^ nr,r,n 

■ , , , , , , r tor a supernova at optical wavelengths (Fian et al. 2005; 

mg to the long burst class, likely originate m the deaths ol 

very massive stars. These massive stars end their lives in a 



Mazzali et al. 2006; SoUerman et al. 2006), while the origin 

^ ^, , 1 • 1 • 1 of the thermal X-ray emission remains an open issue. Shock 

type Ib/c supernova (SN ) explosion, whose signature can be 11, , ■ , 1 , 1 11 

' , , , , , ^tTt-, breaitout must occur m these systems, but its observable 

seen at optical wavelengths alongside the GRB emission m . , . , ,, , 

, signature is not well known, 
a number ol nearby examples with well sampled light curves 

and/or spectroscopy (e.g. Hjorth & Bloom 2011). Supernova Supernova shock break out may have been observed in 

signatures have also been claimed in the high energy emis- a few 'ordinary' supernovae at ultraviolet and optical wave- 

sion of the extremely long, subenergetic and rather unusual lengths (e.g. Schawinski et al. 2008; Gezari et al. 2008,2010; 

gamma-ray burst GRB 060218 in the form of a thermal X- Ofek et al. 2010), but thus far these have been SNe of Type 

ray component present in the spectrum over the first thou- II, and therefore not linked to the GRB systems. Shock 

sand seconds (Campana et al. 2006) that some attribute to break out is the interpretation given to a short burst of X- 

shock breakout of the supernova as it emerges from the star rays observed from a supernova not associated with a GRB, 

(e.g. Waxman, Meszaros & Campana 2007; Nakar & Sari SN2008D, serendipitously caught in the act of exploding 

(Soderberg et al. 2008; Chevalier & Fransson 2008; Modjaz 
et al. 2009; Suzuki & Shigeyama 2010; Balberg & Loeb 2011; 

* Royal Society Dorothy Hodgkin Fellow, mailtoirlcslQle. ac.uk Couch et al. 2011 among others). These serendipitous Swift 
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data on SN2008D provided the earliest ever observations of 
a supernova, seen to rise sharply at X-ray energies followed 
by an exponential decay. In contrast to GRB 060218, the 
X-ray spectrum was not thermal but well described by a 
power law throughout. Shock breakout can lead to either 
thermal or non-thermal spectra depending on the condi- 
tions in the shock and the progenitor star (e.g. Waxman 
et al. 2007; Nakar & Sari 2010; Nakar & Sari 2012), but 
there are other emission sites from which similar spectra 
are expected, namely the cocoon surrounding the jet (e.g. 
Pe'er, Meszaros & Rees 2006). Furthermore, it is unclear 
whether the energies required to produce the observed ther- 
mal X-ray components in GRBs can be obtained through the 
shock breakout process (e.g. Ghisellini, Ghirlanda & Tavec- 
chio 2007). These authors instead invoke the central engine 
to produce the observed X-ray emission. Recently, the no- 
tion of a 'failed GRB' has been suggested for low-luminosity 
GRBs like 060218 (also known as X-ray flashes), in which 
the GRB jet never breaks out of the star yet the choked 
jet powers a relativistic shock breakout which emits in X- 
rays and/or gamma-rays (Bromberg, Nakar & Piran 2011; 
Bromberg et al. 2012). 

A thermal X-ray component has also been noted in 
the prompt to afterglow transition phase of GRB 100316D 
which displayed similar, unusual high energy properties to 
GRB 060218 and is spectroscopically associated with an op- 
tical supernova (Starling et al. 2011). The X-ray spectra of 
these two bursts before about To-|-1000s cannot be fit with 
the absorbed power law or cut-off power law model that de- 
scribes the vast majority of GRBs. The addition of a ther- 
mal, blackbody component with a temperature of '^O.l keV 
significantly improves the fits. 

So it appeared that very long-duration, sub-energetic, 
nearby GRBs with associated supernovae harboured a 
source of thermal X-rays which regular GRBs did not. A 
subsequently reported detection of a thermal X-ray compo- 
nent in GRB 090618 (Page et al. 2011) reopened the debate, 
however, on the origins of such emission. GRB 090618 is a 
'typical' GRB in many respects and yet an absorbed power 
law model was not sufficient to describe its early X-ray spec- 
trum. This GRB has a photometric supernova association 
(Cano et al. 2011a) and lies relatively nearby, at z = 0.54, 
considering the GRB redshift distribution which currently 
stretches from 0.001 < z < 9.4, peaking around z — 2.2 
(e.g. Fynbo et al. 2009; Jakobsson et al. 2012). The large ra- 
dius and necessarily high luminosity of the blackbody seen 
in 090618 is a real challenge to many shock breakout models, 
and in Section [6] we will also consider a model for emission 
from a cocoon surrounding the jet. We note that the early 
X-ray spectrum of GRB 101225A may also be better fit with 
the addition of a thermal component (Campana et al. 2011a; 
Thone et al. 2011), but the redshift is unknown and its clas- 
sification as a GRB has not been confirmed however. 

Whatever its origin, an additional, likely thermal com- 
ponent is required to explain the early X-ray emission of 
at least three GRBs (Campana et al. 2006; Starling et al. 
2011; Page et al. 2011). If it is indeed related to the super- 
nova, and we expect all long GRBs to originate in the deaths 
of massive stars, we should see this component in all long 
GRB which are close enough and the GRB emission is faint 
enough for it to be detected. We note, however, the curious 
cases of long GRBs 060505 and 060614 which possessed all 



the attributes conducive to a SN search, yet no SN could be 
found to very deep limits (e.g. Delia Valle et al. 2006; Fynbo 
et al. 2006; Gal- Yam et al. 2006), highlighting the need for 
a better understanding of the GRB-SN connection. 

We attempt here a systematic search for thermal X-ray 
signatures in a sample of supernova-associated long GRBs 
observed promptly with Swift's X-ray Telescope (XRT, Bur- 
rows et al. 2005; Gehrels et al. 2004). We describe the sample 
in Section [2] and outline the analysis method in Section [S] 
Results for each individual source are presented in Section 
|4] and a supernova-less GRB is discussed in Section [5] In 
Section |S] we summarise the overall findings and outline the 
main caveats. We look at the prompt, afterglow, supernova 
and host galaxy emission of our sample for any relationship 
with the presence/absence of thermal X-ray components and 
we speculate on the origin of GRB-SN thermal X-ray emis- 
sion. In a second, related paper we extend this to all Swift 
GRBs with redshifts, and show simulations intended to re- 
veal the conditions under which the blackbody components 
we are finding can be reliably recovered (Sparre & Starling 
2012, hereafter Paper II). Our conclusions are summarised 
in Section [T] 



2 SAMPLE SELECTION 

We begin our systematic search for additional, thermal X- 
ray components by defining a sample of those GRBs with 

(i) optical SN signatures, either spectroscopic or photo- 
metric, 

(ii) low redshifts, z ^ 1, and 

(iii) sufficient Swift XRT Windowed Timing mode data 
up to 1000 s after the trigger to create a spectrum with 
count rate of a few count s~^ or more. 

The resulting sample includes ten GRBs: 060218, 060729, 
070419A, 080319B, 081007, 090618, 100316D, 100418A, 
101219B and 120422A. We also include GRB 091127, for 
which Swift observations did not start promptly and the 
WT data are beyond 1000 s, since this source shows a very 
clear supernova signature in the optical (SN2009nz, Cobb 
et al. 2010; Berger et al. 2011; Vergani et al. 2011), and 
deserved to be investigated here. It has also been suggested 
that 091127 could be a higher redshift analogue of the suben- 
ergetic bursts such as 980425 and 060218 (Troja et al. 2012). 
Unfortunately we could not apply the same argument to 
GRB050525A, associated with SN2005nc (Delia Valle et 
al. 2006), because the first Swift XRT observations of this 
source from which we could extract a spectrum were both 
late-time and in PC mode due to the faintness of the source 
at that time. For completeness we note two further Swift 
GRBs with SN associations for which prompt XRT spec- 
tra were not collected: 050824 (SoUerman et al. 2007) and 
111211A (de Ugarte Postigo, Thoene & Gorosabel 2012). 
The sample and their redshifts are listed in Table [T] 



3 OBSERVATIONS AND METHOD 

We obtained the Swift XRT reduced spectra for all of our 
GRBs from the UK Swtft Science Data Centre XRT GRB 
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Table 1. Properties of the selected GRB-SN sample. Associated supernovae (P = photometric detection, S = spectroscopic detection to 
be announced), redshift z and Galactic column density TVjj Gal (LAB HI Survey, Kalberla et al. 2005). We also present here the intrinsic 
column density A'^H.int and fit statistic, after fitting the late-time PC mode spectra from the Swift XRT ORB Repository using Cash 
statistics (upper limits are shown in italics). 

*Footnote to table Redshift references: 1-Mirabal & Halpern 2006; SoUerman et al. 2006, 2-Fynbo et al. 2009, 3-Berger et al. 2008, 
4-Cenko et al. 2009, 5-Cucchiara et al. 2009; Thoene et al. 2009, 6-Starling et al. 2011, 7-de Ugarte Postigo et al. 2011a, 8-de Ugarte 
Postigo et al. 2011b; Sparre et al. 2011, 9-Tanvir et al. 2012. 

Supernova references: A-Pian et al. 2006; Sollerman et al. 2006, B-Cano et al. 2011a, C-Hill et al. 2007, D-Bloom et al. 2009, E-Della 
Valle et al. 2008; Soderberg, Berger & Fox 2008, F-Cobb et al. 2010, G-Starling et al. 2011, H-HoUand et al. 2010; de Ugarte Postigo 
et al. in preparation, I-Sparre et al. 2011, J-Wiersema et al. 2012; Melandri et al. 2012; Schulze et al. in preparation. 



ORB 


SN 


iVn.Gai 

xlO^^ 

(cm-2) 


z 


Refs* 


A^H.int 
Xl022 
(cm— 2) 


Cstat (dof) 


060218 


2006aj 


0.094 


0.0331 


1,A 


0.58±0.07 


360 (326) 


060729 


P 


0.045 


0.5428 


2,B 


0.15±0.02 


700 (679) 


070419A 


P 


0.024 


0.9705 


2,C 




6.4 (4) 


080319B 


P 


0.011 


0.9382 


2,D 


0.12±0.06 


493 (537) 


081007 


2008hw 


0.014 


0.5295 


3,E 


0.7±0.2 


193 (276) 


090618 


P 


0.058 


0.54 


4,B 


0.30±0.04 


573 (673) 


091127 


2009nz 


0.028 


0.49 


5,F 


0.11±0.07 


399 (396) 


100316D 


2010bh 


0.070 


0.0591 


6,G 


6+^-2 


25 (21) 


100418A 


S 


0.048 


0.6239 


7,H 


U.<i_Q 2 


32 (35) 


101219B 


2010ma 


0.031 


0.5519 


8,1 




138 (165) 


120422A 


2012bz 


0.037 


0.28 


9, J 


"■'J'^ ' -0.027 


67 (84) 



PL 
PL+BB 
W, fixed 



Late-time well constrained (Table 1)? 
yes no 



PL 
PL+BB 
N,, free 



Addition of BB improves the fit? 
yes report m Table 3 no 



eCPL 
N,, fixed and/or free 



<10 keV break in 
BAT-XRT fit (Table 2)? 



yes 



BKNPL 
N fixed and/or free 



report in Table 4 



report eCPL 
fits in Tabie 4 



Figure 2. Flow chart summarising the sequence of model fits 
applied to the GRB-SN WT mode spectra. 



RepositorjQ. As a first step, we chose to use the complete 
time-averaged WT mode spectrum for all sources (except- 
ing GRB 090618, whose WT data we truncated after the 
steep decay phase according to Page et al. 2011), ensuring 
that the data were early-time. Each spectrum was extracted 
in an identical fashion, using the Swift software version 3.8 
and following the method of Evans et al. (2009). The lat- 
est calibration data files (release 20120209) were adopted. 
The method, including correction for pile-up, hot columns, 
vignetting and exposure maps, is described fully in Evans et 
al. (2007,2009). 

Many GRB time-averaged WT spectra occur during a 
steep decay phase in the fight curve and the time-averaged 
WT spectrum should be sufficient to determine the need for 
a thermal component, but some span more complex tem- 
poral behaviour. We examined the XRT GRB Repository^ 
fight curves for the GRB-SN sample to identify any large 
flares or other spectral anomalies to be avoided and to 
choose, where necessary, intervals over which to extract 
time-sliced spectra. This provided the time ranges reported 
in Table 2, column 2. For those spectra which spanned multi- 
ple light curve decay phases (and therefore temporal slopes) 
we divided the spectrum into time slices where possible (Fig- 
ure [T]). Time slices were chosen to include only one 'decay 
phase' in the temporal behaviour to avoid coadding seg- 
ments during which the major contributing emission com- 
ponents change. 

All WT spectra were grouped such that a minimum of 
20 counts occupied each bin. This allowed statistics to 
be used. 

We also obtained the late-time (>12h) Photon Counting 
(PC) mode spectrum for each GRB, and fit these with an 
absorbed power law model typical of a GRB afterglow, to 



|ww« ■ swift ■ ac . uk/xrt_live_cat| 
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Figure 1. Swift XRT observed count rate ligiit curves of the GRB-SN sample, all shown on the same scale and adapted from the XRT 
GRB Repository. When the WT data were time-sliced, we have indicated the dividing times with grey vertical lines. A colour version 
appears on-line: in cyan are WT settling mode data, blue are WT mode data (during which we fit to the spectra) and red are PC mode 
data. 



obtain a value for the intrinsic X-ray column density, A'h,!!!!, 
which we then fixed in our fits to the WT spectra (listed 
in Table [TJ . Note that this makes the usual assumption of 
constant column density throughout our observations (e.g. 
no detectable ionisation of the circumburst medium) . While 
this is true in general, there are exceptions for which variable 
columns have been (sometimes tentatively) proposed (e.g. 
GRB 011121 - Piro et al. 2005 and GRB 050904 - Boer et al. 
2006; Campana et al. 2007; Gendre et al. 2007; Watson et al. 
2007, but see also Butler & Kocevski 2007). AU of our late- 
time PC mode fits were carried out on essentially ungrouped 
spectra (minimum of 1 count per bin) using Cash statistics 
(Cash 1979) given the low S/N of a minority of the spectra. 
Prompt GRB emission as measured in the hard X- 



rays/soft gamma-rays is usually well described either by a 
power law, or by a cut-off power law or Band function when 
the spectral peak lies within the measured energy range. 
Early X-ray spectra may fall in a transition period between 
the prompt and afterglow regimes, leaving the possibility 
that the spectral peak may enter the XRT band. This is es- 
pecially relevant for the low-luminosity sources which have 
lower energy prompt emission. We performed joint fits to 
both XRT and Burst Alert Telescope (BAT) 15-150 keV 
spectra where both were available, to more robustly lo- 
cate the spectral peak energy or other spectral breaks. The 
BAT data were extracted using the standard Swift reduc- 
tion pipelines and the joint spectra were fitted with ex- 
ponentially cut-off power law models absorbed by a fixed 
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Galactic+intrinsic column (Table [T]), with and without a 
blackbody. Results of these fits are reported in Table (2] 

All spectral analysis was carried out using the X-ray 
spectral fitting package Xspec version 12.7.0 (Arnaud 1996). 

The sequence of model fits to the WT spectra (both 
time-averaged and time-sliced) is summarised in Figure [2] 
To test for the presence of a thermal component we per- 
formed several fits to the early-time WT spectra of each 
of our GRB-SN sample sources, beginning with a power law 
(PL), then adding in a blackbody (BB) at the source redshift 
to represent a thermal component and testing for improve- 
ment using the F-test (where F-statistic < 0.001 indicates 
a significant improvement in fit). Where a BB improved the 
fit, we continued by removing the BB and replacing the 
PL with an exponentially cut-off power law (eCPL). If the 
BAT-XRT joint fits (or literature on the prompt gamma-ray 
data alone) indicated that a spectral peak may be present in 
the XRT band (Table [2]), we also tried a broken power law 
(BKNPL) model. These models mimic the spectral shape 
when a break is present, for example the spectral peak en- 
ergy i5pk as often seen in the prompt emission gamma-ray 
spectrum and shown to be present in the underlying X-ray 
continua of, among our sample, GRBs 060218 (e.g. Kaneko 
et al. 2007) and 100316D (Starhng et al. 2011). 

Galactic absorption, A'^e.Gai, was kept fixed at the 
value from the Leiden Argentine Bonn (LAB) HI Survey 
(Kalberla et al. 2005). Additional absorption at the source, 
A'n.int, was included with the Xspec phabs model using the 
optically-measured redshifts given in Table [1] Solar abun- 
dances according to Wilms, Allen & McCray (2000) and 
cross-sections derived by Verner et al. (1996). The intrinsic 
absorption was initially fixed at the late-time-derived value 
(Table [T]), and allowed to vary in all cases where the 90% 
confidence error on the late-time A'n.int measurement was 
^ 2x the value and/or there were fewer than 100 degrees 
of freedom (dof). A grid of fixed A^H.int values was used 
for GRB 100316D because the late-time PC data were too 
sparse to obtain an accurate spectral fit. 

To further assess the requirement for the newly reported 
thermal components we performed monte carlo simulations 
of an absorbed PL with fixed TVh. We created 10000 simu- 
lated spectra per GRB using this model, fitted these with 
both power law and BB-I-PL models and deduced the per- 
centage of these which recover a (non-existent) thermal com- 
ponent by chance. For this number of trials we can at best 
set a lower limit on the significance of a BB feature of > 4(t, 
allowing low significances (e.g. < 3a) to be identified for 
a modest amount of compute time. Details of the monte 
carlo technique can be found in Hurkett et al. (2008) and 
our results are given in discussions of the individual bursts 
(Section 131). 



4 RESULTS 

The X-ray fitting results for each source are given hereafter 
and in Tables [3] and |3J and Figure [3] We also refer here to 
BAT-XRT fit results, listed in Table [1 GRBs are fisted in 
date order per subsection. 

We find that additional soft X-ray components, de- 
scribed by BBs, are required and favoured in 4 of the GRB- 
SN sample we investigate, with four further possible BBs. 



Table 2. BAT-XRT fits listing the best-fitting model from eCPL 
or eCPL+BB and fitted peak energy, Epj^. For GRB 100316D we 
provide estimates for Epy^ using A''H,int = 0, 0.6, 0.91 and 1.8x10^^ 
cm"'^ in that order. 





'DAT' ■VDT' 

ovcrlSyp {s^ 


model 


^ ,2 / £ 

X /aoj 


TP 

(kcV) 


060218 


160-297 


eCPL-l-BB 


289/296 


37+35 

-13 


060729 


130-347 


eCPL 


268/233 


unconstrained 


060729 


130-160 


eCPL 


152/140 


1 1+°-^ 


060729 


195-347 


too weak to 


fit 




070419A 


119-145 


eCPL 


95/113 


unconstrained 


080319B 


123-182 


eCPL 


267/270 


7+8 
'-4 


081007 


no overlap 








090618 


125-233 


eCPL-fBB 


208/195 


7+11 
'-4 


091127 


no overlap 








100316D 


144-240 


eCPL-l-BB 


234/183 


21±3 


100316D 


144-240 


eCPL-l-BB 


176/183 


19t^ 


100316D 


144-240 


eCPL-l-BB 


173/183 


19+^3 


100316D 


144-240 


eCPL+BB 


189/183 


21+27 
^-^-11 


100418A 


no overlap 








101219B 


153-540 


too weak to 


fit 




120422A 


no overlap 









Three of the sources which we find require a BB are con- 
firmations of previous findings: to avoid unnecessary repeti- 
tion we do not treat these sources in detail but simply show 
that our findings are consistent. The remaining BB exam- 
ple, GRB101219B/SN2010ma, is shown in detail here for 
the first time. We caution that the thermal emission lumi- 
nosities and radii we derive here are only valid if the single 
temperature blackbody is a reasonable approximation and 
the gas is in thermal equilibrium. 

4.1 GRB-SNe with thermal X-ray emission 
GRB 060218 

We recover the thermal X-ray component reported in 
the literature for GRB 060218 (e.g. Campana et al. 2006) 
with high confidence (F-test statistic = 10"^"^). The 
BB-like emission contributes 13% of the X-ray fiux in the 
time-averaged WT spectrum. An eCPL cannot replace 
the addition of a thermal component, and in fact a break 
in the underlying power law in addition to the BB is 
required. This confirms previous work, where it is also 
seen in time-sliced analyses (e.g. Campana et al. 2006; 
Kaneko et al. 2007). The absorbed BKNPL+BB spectrum 
we measure has properties consistent with those published 
elsewhere, validating the method used here. The 0.3-10 
keV luminosity we derive for this BB of ^ 10*^ erg s~^ 
is, together with that of GRB 100316D, the lowest in the 
sample, lying approximately 2 orders of magnitude below 
all the other BB luminosities (see Tables [3] and O but note 
the differing time ranges over which this is calculated). 

GRB 090618 

Our PL and PL-fBB fits for GRB 090618 are consistent 
with the findings of Page et al. (2011) showing an evolving 
BB component, with a restframe starting temperature of 
0.9 keV cooling to 0.5 keV at later times, and an underlying 
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spectral softening seen via the power law photon index. 
The BB is significant at the 99.8% confidence level in the 
time-averaged spectrum, and at 99.998% in the post-flare 
spectrum. Our BAT-XRT fits indicate a spectral break ^ 3 
keV in the overlapping times of 125-233 s since trigger. 
Page et al. (2011) found that this could be constrained 
further when adding in Fermi GBM data, to lie outside the 
XRT energy band at 13±3 keV. We do not perform further 
time-sliced analyses since this exists in previously published 
works. For the purposes of comparison, however, we note 
that we have fixed A'^e.int at a slightly different value (using 
the late-time XRT spectral fit) than did Page et al. (who 
used the XMM-Newton spectrum, see also Campana et al. 
2011b). 

GRB 100316D 

This burst is the subject of detailed analyses in Starling et 
al. (2011) and Olivares et al. (2012) and so only the time- 
averaged early spectrum is treated here. The intrinsic X-ray 
column density cannot be adequately constrained from the 
late-time spectrum, which returns A'n.int = 0.6lJ;g x 10^^ 
cm~'^, so we adopted 4 values for A'n.int of 0, 0.6, 0.91 and 
1.8 xlO^^ cm"^ (from Table[T]or Starling et al. 2011) to fix 
in turn. Whichever value we adopt, excepting A''H,int = 0, 
a BB is required by the F-test at 99.6% confidence. We 
also performed fits with A^n.int free to vary and in this case 
inclusion of the BB improves the fit at the 98.7% confidence 
level. A zero intrinsic column density does not provide an 
acceptable fit with either PL or PL-I-BB models, hence we 
conclude that A^H.int = is excluded. 

The spectral peak energy lies outside the XRT band 
until approximately 240 s in our BAT-XRT fits, and was 
shown to remain just outside the XRT band up until at 
least 737 s by Starling et al. (2011, their Figure 4 and 
Table 3). The BAT-XRT fits also all required a BB. All our 
results hint that this source does indeed have a thermal 
component in its early evolution, additionally supported by 
circumstantial evidence from its remarkable similarities to 
060218, and evidence for thermal emission moving into the 
optical bands following the expected cooling (Olivares et al. 
2012). Because of the difficulty in assigning a correct A''H,int 
value we do not continue with further more complex fits to 
this source. The important role that the intrinsic column 
density plays in determination of any BB component is 
discussed in detail in Paper II. 

GRB 101219B 

The time-averaged spectrum cannot be modelled with a sin- 
gle absorbed PL, and inclusion of a BB with temperature 
kT = 0.2 keV and luminosity ^10*^ erg s~^ improves the 
fit very significantly. Monte carlo simulations (as described 
for GRB 100316D in Starling et al. 2011 and more generally 
in Hurkett et al. 2008) show that the BB component in this 
spectrum is > 4(7 significant. We tested an absorbed eCPL 
model as well, and this provided no improvement over the 
PL model. 

Dividing this time interval into three we studied the 
time-sliced spectra in the same manner. From 160-256 s 
since trigger the BB improves the fit at 99.99% confidence, at 
256-320 s the confidence level is 99.79% and over 320-544 s 
there is improvement at very high confidence, >99.99%, as 
found in the time-averaged spectrum. The BB temperature 



decreases with time over the full duration of the WT ob- 
servations, from 0.30±0.06 to 0.18±0.02 keV, while the lu- 
minosity approximately doubles and its contribution to the 
unabsorbed 0.3-10 keV observer-frame flux goes from 6% to 
25%. At the same time the power law photon index steep- 
ens, seen in a number of GRBs during sharp decays in flux 
in the prompt to afterglow transition (Zhang et al. 2007). 

All the WT data were taken during an overall steep de- 
cay phase in the light curve. This decay continues after the 
switch to PC mode, and it is therefore interesting to extract 
the flrst orbit PC spectrum covering 544-1056 s and flt the 
same model. Statistics are poor, and we group the spectrum 
such that a minimum of 15 counts lie in each bin. If we per- 
form a PL-I-BB flt we obtain the following BB parameters: 
kT = 0.10±0.02 keV, F = 1.6±0.1, X^/dof = 16.2/16 show- 
ing a continuation of the cooling measured in earlier data. 
The BB by this time seems to be fading again, with luminos- 
ity similar to that measured in the second time-sliced WT 
spectrum (but with larger error bar), now contributing 34% 
of the flux. In Figure |4] we show the BB-I-PL spectral flts 
to each of the time-sliced spectra described above, and in 
Figure [5] the evolutionary trends of key spectral parameters 
over 160-1056 s are shown. 

4.2 Possible thermal emission candidates 
GRB 060729 

A BB+PL is a signiflcantly better flt over a PL for the 
time-averaged WT spectrum. However, we note a flare in 
the light curve, which typically produces a hardening in 
GRB spectra, so we split the spectrum in two segments, 
pre- and post-flare. The pre-flare spectral flt does not im- 
prove with the introduction of a BB component, though it 
is also not well flt with the absorbed PL model, while the 
post-flare spectrum does favour such an addition. It is clear 
that the time-averaged spectrum follows most closely the 
post-flare spectrum which contains more counts than pre- 
flare, and we note that the hardness ratio (HR) on the XRT 
GRB Repository^ shows much spectral softening during the 
pre-flare spectrum. Monte carlo simulations show that the 
BB component in the post-flare spectrum is > 4(t signifl- 
cant. Grupe et al. (2007) performed a detailed analysis of 
GRB 060729. They also found that an additional BB com- 
ponent improves the flts for the later WT data: although 
they consider very different time-slices from us, they state 
that a BB is required after 150 s, which is consistent with 
the presence of this component in our second time interval of 
195-352 s. Their range of BB temperatures is also consistent 
with the mean value determined here. 

Over 130-160 s, during the pre-flare spectrum, we have 
overlap between the BAT and XRT coverage. Fitting the 
joint BAT-XRT spectrum reveals a spectral break around 1 
keV. Accounting for this using the broken power law model 
we flnd the pre-flare XRT spectrum can be well fltted with a 
break that is consistent with the BAT-XRT joint flt result. 
Post-flare we have too little BAT data to perform joint BAT- 
XRT spectral flts. Fitting the XRT post-flare spectrum with 
a BKNPL instead of a PL-l-BB we flnd that while it remains 
an improvement over a single PL model, the flt statistic is 
not as good as the addition of a BB and although _Ebk moves 
to lower energies as expected, the spectral slopes are very 
different from those found pre-flare. 
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It is perhaps unusual, given previous examples (al- 
though few in number), that a BB is not apparent in the 
earliest XRT data, but is measureable in subsequent spec- 
tra. This is a bright afterglow, which may drown out low 
luminosity components; the BB is detectable in GRB 090618 
which lies at the same redshift, has a similar shaped light 
curve and is even brighter in the first 200 seconds, but the 
BB luminosity we find for 060729 post-flare is only 1/3 of 
that in 090618. For 060729 we can only conclude that there 
is significant curvature in the spectrum which may be a 
hidden BB, but its detection would be complicated by the 
flaring episode and presence of the spectral peak cascading 
through the energy band and we are unable to draw firm 
conclusions. 

GRB 081007 

A model comprising a PL+BB is a better fit than a single 
PL for this source with significance 99.99%. In monte 
carlo simulations of an absorbed power law, we found a 
chance improvement in the fit with a BB in 6/9990 trials, 
reducing the significance of the result for 081007 in Table [3] 
to 99.94%. Our eCPL fits, with and without a BB, did not 
improve the fit statistic because the fitted cut-off energy 
was unconstrained, tending to far beyond the XRT energy 
band consistent with the single PL fit to prompt Swift BAT 
and Fermi GBM spectra (Markwardt et al. 2008; Bissaldi, 
McBreen & Connaughton 2008). The XRT spectrum has 
relatively few total counts, resulting in substantial errors 
on the power law slope when a BB is included in the fixed 
A'H.int fit (Figure We note that A^n.int is reasonably 
large, which may present difficulties in identifying any 
BB component (see Paper II). While we only show the 
fixed A'^H.int fits in Table O (as the late-time-derived A''H,int 
is known to within ~30%), allowing A'n.int to vary also 
resulted in an acceptable fit with a steeper power law 
photon index and no need for a BB. Intrinsic column 
density in this case increases, but remains just consistent at 
the 90% level with the late-time fit. For this reason, we can 
only report here a possible, unconfirmed BB component. 
This is one of the fainter sources in WT mode in our sample. 

GRB 100418A 

The addition of a BB when the column density is fixed 
resulted in a very significant improvement in the fit 
statistic and monte carlo simulations show that the BB 
component in this spectrum is > 4(t significant. However, 
the uncertainties on A^n.int are large and allowing this 
to be a free parameter in the PL model also provided 
a significant improvement to the fit while increasing the 
column density from the fixed value of 2x10^^ to (6±l)x 
lO'^^ cm~'^. The addition of a BB at that stage provided 
another small improvement in the fit, but interestingly the 
column density then reverts to around its originally fixed 
value. An eCPL model with fixed A'n.int is not an adequate 
fit. With A'H,mt left free (again converging to ~6xl0^^ 
cm~^) this improves, but the cut-off lies way outside the 
considered energy band meaning we are essentially just 
fitting a power law as before. This is consistent with the 
single PL fits reported for prompt BAT data (Ukwatta et 
al. 2010). Trying instead a BKNPL, we find an adequate fit 
with Ebk ~ 0.6 ± 0.1 keV. It is very difficult, however, to 
constrain a break at such low energy. If we then allow A''H,iiit 



to vary it is apparent that A^n.int and Fi are degenerate 
(Figure [3]) so we do not consider this result further. In 
conclusion, the PL+BB model with or without A'n.int 
free has the best fit statistic and simple alteration of the 
PL model with free A^H.int parameter values alone does 
not seem to achieve this but also results in an acceptable fit. 

GRB 120422A 

The WT spectrum cannot be modelled with a single ab- 
sorbed PL. Inclusion of a BB significantly improves the fit 
and monte carlo simulations show that the BB component 
in this spectrum is > 4(t significant. However, residuals still 
remain and Xrcducod ~ If the intrinsic absorbing column 
is allowed to vary, as shown in Table [3] we still prefer the 
inclusion of a BB both due to improvement in the fit and to 
alleviate the need for a very steep PL photon index. The fit 
statistic for the BB-I-PL model with An.int free, however, is 
of the same order as that when TVeant is fixed, and A^n.int is 
poorly constrained. The restframe BB temperature is of or- 
der 0.2 keV with a luminosity of ~10''^ erg s~^, contributing 
50-60% of the unabsorbed 0.3-10 keV fiux. 

We do not have any BAT-XRT overlap or detection 
in the literature of a spectral peak energy in the prompt 
phase. An eCPL model does not provide an improvement 
over the single PL, with the peak energy coverging to a 
value high above the XRT band. Instead we tried a BKNPL 
- the fit statistic is acceptable, though not quite as good as 
the BB-I-PL fits. It resulted in a break energy of 0.6 keV and 
a poorly constrained Fi. The BB+PL model is preferred for 
this source, but a BKNPL and/or increase in the very low 
intrinsic column cannot be ruled out. We note that, using a 
cut-off PL model, Zhang et al. (2012) suggest this source is 
inconsistent with shock breakout models and suggest instead 
an engine-driven GRB with steep X-ray decay caused by the 
curvature effect. 

4.3 GRB-SNe with no detectable thermal X-ray 
emission 

GRB 070419A 

In addition to the time-averaged WT spectrum we created 
a further time-selected spectrum for GRB070419A to avoid 
early fiaring. In both spectra a fixed A'^e.int PL-fBB is a 
better fit than a single PL with F-test statistic of ~5e-3. 
The BB temperatures of 0.1-0.2 keV we measure are similar 
to those found in GRB 060218, while the BB luminosity is 
1000 times greater. However, the intrinsic column density is 
not known to any reasonable degree of accuracy due to poor 
S/N in the late-time spectrum so we also performed a fit 
with A^H.int allowed to vary. In these latter fits we find that 
an absorbed PL is a perfectly adequate description of the 
data: the fit statistic improves significantly and the intrinsic 
column density converges to a lower value than previously 
fixed at while still remaining consistent with AfH.iatotimc (as 
expected given its large error bars). 

BAT-XRT fits in the temporal overlap of 119-145 s 
(covering some of the time-averaged spectrum but only 
during the flaring episode) did not constrain any spectral 
break. An eCPL fit to the XRT spectra after the fiare 
shows that with fixed A^n.int we can obtain a very good fit 
with peak energy i?pk around 4-7 keV (no BB needed) and 
Fi steepens a little compared with the PL-only fits, while 
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when A^H.int is left free to vary i?pk cannot be constrained 
since there are hkely too few photons in the energy bins 
above a few keV. In either case the eCPL fits are of similar 
fit statistic to the PL+BB in the fixed A''H,int case or the 
PL in the free A^'h int casG, thus we conclude that there is 
no evidence to support the presence of a thermal X-ray 
component in this source. This is the most distant GRB-SN 
in our sample at z ~ 0.97. 

GRB 080319B 

There are a lot of WT data for this bright burst, therefore 
we created 3 time-sliced spectra (split by light curve decay 
segments, Fig[T]) in addition to the time-averaged spectrum 
and fit all four. We conclude here that no BB is required. 
However, the single absorbed PL model is a poor fit to the 
data from the time-averaged and final time-slice spectra, 
which resemble each other closely since the latest time-slice 
we have chosen contains a much larger portion of the total 
WT mode exposure time. We therefore tried fitting all of 
our 080319B spectra with the eCPL model. This marginally 
improves the fits overall, but greatly improves the fit to the 
time-averaged and final time-slice spectra, with the peak 
energy seeming to cascade downwards in energy from about 
7 to 4 keV. For the final time slice we also tested a BKNPL 
because the peak energy fell in the mid range of our spectral 
energy band. This further improved the fit and brought the 
spectral break to 1-2 keV. In all these cases we checked the 
result of letting A^H.int vary, and the value did not stray far 
from its late-time measurement. In summary, the early X- 
ray data of 080319B can be fit with a set of power laws 
with a spectral break that moves through the X-ray band 
to lower energies becoming clearly measurable between 704 
and 1744 s. Our BAT-XRT fits support the presence of a 
low energy spectral peak, measuring 3 ^ iJpk ^ 15 keV at 
123-182 s. 

Racusin et al. (2008) present a two-jet model to ex- 
plain the shape of the X-ray afterglow temporal decay. We 
note that the portion of the light curve we sample here is 
dominated by the narrow jet in their model (the wide jet 
contributes approximately 1/100 of the flux in their Supple- 
mentary Figure 7) and does not undergo a temporal break 
until after our observations. In our search for a thermal com- 
ponent we should therefore not be confused by temporal 
breaks or interplay between the proposed two jet compo- 
nents according to their model. These authors require the 
cooling break for the narrow jet forward shock to lie above 
the X-ray band throughout the time covered by our X-ray 
observations, ruling out Vc as the break in our fits. They also 
state that the spectral peak of the narrow jet must cross the 
X-ray band prior to 60 s when XRT observations began, also 
ruling out the break suggested in our fits. 

Bloom et al. (2009), who proposed the existence of an 
optical supernova bump in the light curve (and hence this 
sources' inclusion in this paper), could not fit the early X- 
ray emission from 080319B within the standard model. Their 
preferred model while different from that of Racusin et al. 
also requires that the cooling frequency lies above the X-ray 
regime for the duration of the observable X-ray afterglow, 
again ruling out Vc as the break suggested from our fitting. 
However, these authors focussed on the optical, infrared and 
PC-mode X-ray data and did not attempt to explain the 
early X-ray data, suggesting it is umelated to the early op- 



Figure 5. Evolution of spectral parameters when fitting a 
PL-I-BB model to GRB 101219B time-sliced spectra (Tables l3l5ll. 
The top panel shows steepening of the underlying power law in- 
dex, the second panel from the top shows cooling of the fitted BB. 
The lower two panels show the evolutionary trends of luminosity 
and radius for the BB with errors taken from the BB normalisa- 
tion and temperature (redshift uncertainty is not included) . Time 
since BAT trigger is in the observer frame and plotted on a log 
scale. 
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tical emission since the two have very different behaviours, 
or at least there is a spectral break between those bands 
which moves to higher frequencies with time. We do not see 
a break moving in the direction of higher frequencies, but if 
breaks within the X-ray band, thermal components or other 
spectral curvature were not taken into account then this may 
explain the difficulty in fitting these data. 

We also point out that the redshift of this GRB is 
the second highest included, at 2: ~ 0.94, and that the 
photometric identification of a supernova in this source is 
highly uncertain. 

GRB 091127 

The spectrum of GRB 091127 is adequately fit with an ab- 
sorbed PL. Given the post-3000 s observation start, much 
later than any of the other sample objects, it may be too 
late to see any detectable thermal component. 



5 SUPERNOVA-LESS GRBS 

An intriguing individual case which fits all but the first of our 
sample selection requirements in Section [2] (the requirement 
for an associated SN) is GRB 060614. This source was as- 
signed to the long GRB class (e.g. Xu et al. 2009), while this 



© RAS, MNRAS 000, [THIS] 



A search for thermal X-ray signatures in GRBs I 9 



Table 3. Table of results from absorbed PL and absorbed PL+BB fits to the GRB-SN sample. Column 1 also lists time since BAT 
trigger covered by the spectrum, f denotes a fixed column density. refers to the fraction of the total observer frame 0.3-10 keV ratio 

of unabsorbed flux contributed by the BB component. The total luminosity of the BB component is calculated from the BB normalisation 
and luminosity distance, as given in Page et al. (2011). F-statistic gives the F-Test probability of a chance improvement in the fit when 
adding a BB component (or freeing A^H.int if given in brackets). 



GRBitime 

(s) 


model 


Xl022 
(cm~2) 


(keV) 


F 


^BB 

fx 

(%) 


^BB 
Xl047 
(erg s^i) 


do/ 


Fstat 


060218:159-2783 
060218:159-2783 


PL 

PL-fBB 


0.58 f 
0.58 f 


0.181+0.004 


1.820+0.007 
1.63+0.01 


13 


0.024+0.001 


2826/788 
1528/786 


10-105 


060729:128-352 
060729:128-352 


PL 

PL-I-BB 


0.15 f 
0.15 f 


302+°-"2^ 


2.82+0.03 
2.74+0.07 


20 


16+^ 


272/169 
184/167 


10-15 


060729:128-160 
060729:128-160 
060729:195-352 
060729:195-352 


PL 

PL-fBB 
PL 

PL-fBB 


0.15 f 
0.15 f 
0.15 f 
0.15 f 


52+'''05 
^■°^-0.04 

0.21+0.01 


1.98+0.06 

1 7+0.2 
^■'-0.3 

3.25+0.05 
2.9+0.2 


37 


22+5 


230/69 
230/67 
230/118 
145/116 


10-12 


070419A:112-304 
070419A:112-304 
070419A:112-304 
070419A:112-304 


PL 

PL+BB 
PL 

PL+BB 


1 f 
1 f 

n O1+0.09 

^•^^-0.08 

0.76to:i? 


0.19+0.08 
4+°-2 


9 q+0.05 

^■•^-0.049 
9 ig+0.08 
^•^^-0.10 

2.19+0.07 
2 1+0-1 


11 
6 


^°-50 

31+31 


179/166 
168/164 
167/165 
166/163 


6x10-3 

(7x10-4) 

0.61 


070419A: 150-304 
070419A: 150-304 
070419A: 150-304 
070419A: 150-304 


PL 

PL+BB 
PL 

PL+BB 


1 f 
1 f 

"^-0.09 
Q 7+0.3 

'-0.2 


'^■^°-0.05 


2.42+0.06 
2.31+0.09 

9 90+0. 09 
^■^"-0.08 
9 99+0.42 
^■^^-0.01 


16 


-1 77+490 
^' '-102 


119/111 
108/109 
105/110 
105/108 


5x10-3 
(2x10-4) 


080319B:60-1744 
080319B:60-1744 


PL 

PL+BB 


0.12 f 
0.12 f 


0.93+0.05 


1.67+0.008 
1.67+0.015 






1213/592 
1214/590 





080319B:60-304 

080319B:60-304 

080319B:304-704 

080319B:304-704 

080319B:704-1744 

080319B:704-1744 



PL 

PL+BB 
PL 

PL+BB 
PL 

PL+BB 



0.12 f 
0.12 f 
0.12 f 
0.12 f 
0.12 f 
0.12 f 



1.0+0.4 
0.8+0.1 
0.92+0.04 



1.71+0.03 

1.72+0.04 

1.64+0.02 

1.62+0.03 

1.660+0.008 

1.66+0.02 



420 



+316 
255 



295/266 

288/264 

486/386 

486/384 

1274/572 

1274/570 



0.04 



081007:105-175 
081007:105-175 


PL 

PL+BB 


0.7 f 
0.7 f 


0.31+0.04 


2.55+0.09 
2.1+0.3 


32 


61^3 


45.6/38 
29.7/36 


4x10-" 


090618:125-250 
090618:125-250 


PL 

PL+BB 


0.3 f 
0.3 f 


Q+O-3 


1.45+0.05 
1.42+0.07 


10 


18111 


145/133 
132/131 


2x10-3 


090618:160-250 
090618:160-250 


PL 

PL+BB 


0.3 f 
0.3 f 


51+0'08 
^■°l-0.05 


1.9+0.04 
1.83+0.06 


11 
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113/129 


2x10-5 


091127:3232-3840 
091127:3232-3840 


PL 

PL+BB 


0.11 f 
0.11 f 


6+"-3 


1.77+0.04 
1.75+0.06 






158/176 
158/174 





1003160:144-736 


PL 


f 




0.42+0.02 






3628/491 




1003160:144-736 


PL+BB 


f 


001+°-''3'^ 


0.42+0.02 


<1 




3628/489 




1003160:144-736 


PL 


0.6 f 




1.11+0.02 






826/491 




1003160:144-736 


PL+BB 


0.6 f 


0.72+0.02 


0.80+0.05 


20 


0.046+0.006 


590/489 


10-36 


1003160:144-736 


PL 


0.91 f 




1.33+0.02 






536/491 




1003160:144-736 


PL+BB 


0.91 f 


n. c:o + 0.12 

U-aS-o.oe 


1.24+0.04 


5 


012+0-007 


524/489 


4x10-3 


1003160:144-736 


PL 


1.8 f 




1.77+0.03 






832/491 




1003160:144-736 


PL+BB 


1.8 f 


0.176+0.009 


1.50+0.03 


33 


16+0-05 
'^•^"-0.03 


512/489 


10-52 


1003160:144-736 


PL 


1.07+0.05 




1.43+0.04 






508/490 




1003160:144-736 


PL+BB 


-I 97+0.14 
-^•^'-0.10 


23+0 04 

"•^'-'-0.03 


1 4+0.03 
-^•^-0.04 


8 


026+0-036 


499/488 


0.013 


100418A:64-160 


PL 


0.2 f 




3.27+0.08 






94.9/46 




100418A:64-160 


PL+BB 


0.2 f 


0.23+0.01 


2.7+0.2 


57 


16+4 


34.7/44 


10-10 


100418A:64-160 


PL 


0.6+0.1 




4.2+0.3 






46.0/45 


(10-8) 


100418A:64-160 


PL+BB 


0.281°-^ 


0.22+0.04 


2.9+0.7 


52 


17+12 


34.2/43 


2x10-3 
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Figure 3. On the left hand side we show the unfolded time-averaged WT-mode spectra of GRBs 101219B, 081007, 100418A and 120422A, 
and the post-flare spectrum of 060729, with the best-fitting absorbed BB+PL model with fixed A'^h int overlaid as dotted lines. To the 
right of each spectrum we show the 68/90/99%confidence contour plots for 2 interesting parameters; BB temperature kT, against power 
law photon index F, and kT against intrinsic column density Afjj.int when that was free to vary. 




Energy (keV) Parameter: Pholndex Parameter: nH (1 0=^) 
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Table 3. Table\3\ continued. 



101219B:160-544 


PL 


0.16 f 




1.71±0.05 


292/186 


101219B:160-544 


PL+BB 


0.16 f 


0.20±0.02 


1.34±0.07 14 3.5+05 


159/184 10-25 



101219B:160-256 
101219B:160-256 
101219B:256-320 
101219B:256-320 
101219B:320-544 
101219B:320-544 



PL 0.16 f 

PL+BB 0.16 f 

PL 0.16 f 

PL+BB 0.16 f 

PL 0.16 f 

PL+BB 0.16 f 



0.30±0.06 
0.23±0.06 
0.18±0.02 



1.30±0.07 
l.liO.l 6 

1 37+0.17 q 

2.17±0.08 
1.56±0.12 25 



85.1/76 

2±1 75.1/74 1x10""' 

36.5/42 

2.5±1 26.8/40 2x10-3 
174/89 

4.8±0.7 73.0/87 lO^^'^ 



120422A:128-192 
120422A:128-192 
120422A: 128-192 
120422A:128-192 



PL 0.027 f 

PL+BB 0.027 f 



PL 

PL+BB 



36+° ''^ 
'-'■■'°-0.06 
,+0.14 
-0.11 



0.13-* 



0.204±0.009 2.37 
4.35 



0.17"^ 



2.98±0.06 

+0.27 
0.34 
+0.30 
-0.27 
2 Ql +0.65 
^■yj^-0.62 



65 1.9 



52 2.3 



213/59 

76/57 10-" 

88/58 (10-12) 

+0* 74/56 8x10-3 



Table 4. Exponentially cut-off power law (eCPL) and broken power law (BKNPL) fits to those sources requiring further investigation. 
Column 1 also lists time since BAT trigger covered by the spectrum. 



GRB:time 


model 




ri,2 


-E'pk/bk 


dof 


(s) 




Xl022 (cm-2) 




(kcV) 


060729:128-160 


BKNPL 


0.15 f 


O.Sto l, 2.7±0.1 


-1 17+0.05 
^■^'-0.09 


71/71 


060729:195-352 


BKNPL 


0.15 f 


2.5+°-2, 3.62±0.07 


0.701°°^ 


161/116 


070419A: 112-304 


eCPL 


0.1 f 


^■"^"-0.08 


4.01^4 


166/165 


070419A:112-304 


eCPL 


'^■^'-'-0.09 


2.4±0.1 


4 6+31 


164/164 


070419A: 150-304 


eCPL 


0.1 f 


2.68±0.09 


^■°-0.7 


110/110 


070419A:150-304 


eCPL 


0.8±0.1 


2 4+''-2 
^■^-0.1 


unconstrained 


105/109 


080319B:60-1744 


eCPL 


0.12 f 


1 07+0.02 

-0.03 


3.9±0.2 


713/592 


080319B:60-304 


eCPL 


0.12 f 


1 60+°"^ 

^■""-0.08 


"■"-1.9 


288/265 


080319B:304-704 


eCPL 


0.12 f 


1.43±0.05 




413/386 


080319B:704-1744 


eCPL 


0.12 f 


1.34±0.02 


3.8±0.2 


694/572 


080319B:704-1744 


BKNPL 


0.12 f 


1.41±0.02, 1.87±0.02 


1 40+009 
l-''>J-0.06 


663/571 


080319B:704-1744 


BKNPL 


'^■^'-'-0.03 


l-37to Qg, l-87±0.02 


1 00+0. 07 
1-38-0.05 


661/570 


081007:105-175 


eCPL 


0.7 f 


2.6±0.1 


unconstrained 


45.6/37 


100418A:64-160 


BKNPL 


0.2 f 


3.8±0.2 


0.6±0.1 


40/41 


101219B:160-544 


eCPL 


0.16 f 


1.71±0.05 


unconstrained 


292/185 


120422A: 128-192 


BKNPL 


0.037 f 


0-eto,l, 3.6±0.1 


0.60±0.04 


79/57 



was somewhat controversial given the short spike-|-long tail 
prompt emission profile (e.g. Gehrels et al. 2006). It lies at a 
low redshift, z — 0.125 (e.g. Fynbo et al. 2009), yet no opti- 
cal supernova signature was found to very deep limits (Delia 
Valle et al. 2006; Fynbo et al. 2006; Gal- Yam et al. 2006). 
Mangano et al. (2007) performed detailed fits to the WT 
mode data of 060614, always with a variable intrinsic col- 
umn density. They improved the fit statistic by going from 
an absorbed PL to either PL+BB, eCPL or Band models 
(all equally acceptable), and on inclusion of the BAT data 
during overlap conclude that this was due to the spectral 
peak energy moving through the XRT band from 97-175 s 
after burst. We performed the same fits as in Section [3] on 
the WT time-averaged (97-466 s) spectrum of 060614, for 
consistency and comparison with our GRB-SN sample and 
in the hope that it may shed further light on the nature of 
this unusual GRB. We find that a single absorbed PL pro- 
vides a good fit with no need for further components, when 
fixing the column densities to TVn.Cai = 1.9 x lO^'' cm-^ 
(LAB HI Survey) and iYn.int = 3 x 10^° cm^^ (late-time 
PC mode spectral fit). Allowing the intrinsic column den- 



sity to vary, the PL model is still preferred over a PL+BB. 
These results are in fact consistent with those of Mangano 
et al. (2007): the mean photon arrival time of 250 s for the 
time-averaged spectrum means that iJpk will no longer be 
detectable within the XKT energy band, and a single PL 
will be an adequate fit. 



6 DISCUSSION 

We report a new thermal X-ray component detection for 
GRB 101219B associated with SN2010ma, which we believe 
to be compelling. We also recover the BB components re- 
ported previously in GRBs 060218, 090618 and 100316D 
with consistent parameters. For our remaining 5 BB can- 
didate sources from Table O (F-test statistic ^ 10-^) we 
then omitted times of major fiaring in the light curves, time- 
sliced according to light curve segment, allowed A'^e.int to 
vary when the error from the late-time measurement was 
very large and fit eCPL and BKNPL models. We performed 
joint BAT-XRT fits to all our sample where possible to lo- 
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Figure 4. Early X-ray evolution for GRB 101219B with a PL+BB model fit as listed in Table |3] increasing with time through black- 
red-green. The data are displayed as E'^Fe vs E, equivalent to uF^. The blue data show the addition of the PC-mode spectrum, which 
appears to be a good extrapolation of the earlier evolution. 




Energy (keV) 



Table 5. Measured properties of the final list of BB candidates and possible candidates from the time-averaged spectra as listed in Table 
[3](N.B. Fits to GRBs 100316D, 100418A and 120422A are Afn.int-dependent). For GRB 101219B we also list the time-sliced spectral 
results including PC mode. 

GRB tiogniid Tbb Fbb isB Rbb 





restframe (s) 


(keV) 


(erg cm ^ s ^ ) 


(%) 


xlO**'^ (erg s-l) 


(cm) 


060218 


643.9 


0.181±0.004 


8.10x10"^'' 


13 


0.024 


4x10" 


090618 


114.8 


9+°-^ 


6.530x10-9 


10 


181 


1.5x10^2 


090618 


129.9 


51+0 08 
'^•'^^-0.05 


2.150x10^9 


11 


61 


3x10^2 


100316D 


307.4 


0.72±0.02 


5.00x10""^ 


20 


0.046 


4x10^0 


100316D 




58+''-^2 
U-0»_o.06 


1.30x10-1'' 


5 


0.012 


3xlOio 


100316D 




0.176±0.009 


1.500x10-9 


33 


0.16 


1x10^2 


100316D 




r, r,, + 0.04 


2.58x10-1'' 


8 


0.026 


3x10" 


101219B 


202.4 


0.30±0.06 


6.1x10-" 


6 


1.97 


9x10" 


101219B 


286.2 


0.23±0.06 


6.9x10-" 


9 


2.48 


2x10^2 


101219B 


417.2 


0.18±0.02 


1.13x10-"' 


25 


4.80 


4xl0i2 


101219B 


757.9 


0.10±0.02 


2.6x10-" 


34 


2.50 


7x10^2 



060729 169.8 0.21±0.01 6.31x10-"' 37 22 9x10^2 

081007 88.6 0.31±0.04 2.10x10-1" 32 6 2x10^2 

100418A 62.3 0.23±0.01 3.00x10-1" 57 16 7xl0i2 

100418A 0.22±0.04 3.15x10-1" 52 17 8x101^ 

120422A 122.5 0.204±0.009 3.98x10-1" 65 1.9 3x101^ 

120422A 0.17toQ4 4.49x10-1" 52 2.3 5x101^ 
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cate the spectral peak energy and thereby better diagnose a 
spectral break in the early XRT data. 

We found we could reasonably explain the spectral 
curvature with one of the eCPL or BKNPL models for 
GRB 070419A and the early time-slice of GRB 060729. We 
conclude that four sources may show thermal X-ray emission 
within their early X-ray spectra, but cannot rule out other 
explanations for this curvature. The reasons for this are most 
often because of large uncertainties in absorbing column and 
degeneracy between this and the underlying spectrum, or 
unconfirmed spectral breaks in the spectra. 
The role of X-ray absorbing columns 
The absorbing column along the line-of-sight towards a GRB 
can be substantial, and reduces the observed X-ray flux in 
the region we are interested in for BB detection. It is there- 
fore likely that details in the soft X-ray spectrum will be 
missed if the absorbing column is sufficiently high. The col- 
umn density limit above which this happens is investigated 
in Paper II and must depend on power law index and source 
redshift. We note that the Galactic column towards our 
GRB-SNe has a mean of 4x10^° cm~^, which is insignificant 
when compared with the mean intrinsic measured column of 
~4xlO^^ cm~^ and so should have no effect. Signal to noise 
in many of our spectra is not good enough to keep Ai^H.int 
free and add a BB, hence in the first fits in this work we kept 
A'H,int fixed at the late-time value (Table [1]). This is clearly 
not ideal, and we had to identify the cases with very uncer- 
tain late-time absorbing column densities, Afn.int, and allow 
for this. The BB in 100418A must be deemed uncertain due 
to the errors on A''H,int, while I20422A is consistent with hav- 
ing very little or no intrinsic absorption. For the remaining 
sources, we note that 060218, 081007 and possibly I003I6D 
all have greater intrinsic columns than the mean of the sam- 
ple and, while the BB is clear in 060218 and 100316D, it 
is likely that the high column (A'^H.int = (0.7±0.2) x 10^^ 
cm~^) is affecting our ability to either securely identify or 
rule out a BB component in 081007. 
Spectral evolution 

Spectral evolution at early times in the X-ray production of 
GRBs is common, where the power law slope of the spec- 
trum is seen to steepen through the steep decay phase in 
the light curve (e.g. Zhang, Liang & Zhang 2007). We can 
estimate the effect of spectral evolution on our sources by 
examination of their hardness ratios as created in the light 
curves section of the UKSSDC Swift XRT GRB Repository \ 
and in some cases we split the spectra according to the dif- 
ferent light curve segments (different decay slopes). Harden- 
ing during flares is a particular problem and we have done 
our best to eliminate flaring episodes from our spectra. The 
time-sliced analysis of GRB 101219B shows a softening of 
the power law, but this did not render the BB undetectable 
in the time-averaged data so we expect this effect will not 
be a significant issue. 

The presence of the spectral peak energy in the X-ray 
band during our observations would invalidate our assump- 
tion of a single underlying absorbed power law continuum, 
requiring a spectral break or added curvature for the case of 
a long-duration exposure in which the break energy is seen to 
move. We know that Spk is present in the X-ray spectrum we 
use here for GRB 060218 yet we can still conclusively recover 
the BB component (though it is evident from the poor fit 
statistic that the continuum model remains incorrect, likely 



due to significant spectral evolution). This was the most 
prominent thermal component ever detected, comprising up 
to half of the observed X-ray flux. The exponentially cut- 
off power law fits to the spectra (Table |4]), together with 
information from either prompt gamma-ray fits or our own 
BAT-XRT fits (Table [5]) have allowed us to reasonably as- 
sess the need for a spectral break in each spectrum. Indeed, 
we find that a break is likely to be present in the early X-ray 
data of 070419A, 060729 pre-flare and 080319B (in the lat- 
ter we note that a two-jet model suggested in Racusin et al. 
2008 rather implies a superposition of power laws, but other 
interpretations, namely that of Bloom et al. 2009, require 
only one jet). 
Instrumental effects 

We note that the quasar 3C 273 is, like some of our GRB-SN, 
well described by an absorbed power law plus blackbody and 
may be used to assess the effect of calibration uncertainties 
on our sample. This source is a regular Swift XRT calibration 
target, and was observed in 2005 simultaneously with XMM- 
Newton and RXTE for cross-calibration purposes. The BB 
in 3C 273 has a temperature of kT = 0.1 keV and contributes 
4% of the XRT X-ray flux. The spectral flts for XRT WT 
mode agree to within the expected uncertainties with the 
XMM EPIC and RXTE spectra (we note that for high sta- 
tistical spectra a systematic error of 5C3% was found to be 
sufficient, see e.g. SWIFT-XRT-CALDB-09_vll and subse- 
quent calibration release notes). This gives us confidence 
that the BB components required in the GRB-SN sample 
are real spectral features, and cannot be an instrumental 
effect. 

6.1 Comparison with GRB, SN and host galaxy 
properties 

Although we have small number statistics here with only 
11 sources in our sample, we can examine the GRB, SN 
and host galaxy properties and look for commonalities be- 
tween the source and BB properties. We split our GRBs 
into three 'detection categories': detected BB, possible BB 
and no BB, as they are set out in Section U No correlations 
arise between any of the properties we looked at: z, Tgo, 

•^BATIS-ISO keV, E\so, fobs, start, tniid,rcst, «X , A^H,Gal, A^H,int, 
r, fcTBB, -Fbb/-Fx, iBB, -RbB, F-Statistic, Mv,SN-pcak, SN 

detection technique (spectroscopic or photometric), AfB.host 
when host is detected and BB required (yes, possibly or no). 

There is, unsurprisingly, a general trend with redshift: 
the two most distant sample sources at z ~ 1 do not have 
detectable BB signatures, while BBs in the two very nearby 
sources at 2 < 0.1 are clearly detected. A 0.2 keV BB will 
peak at ~0.1keV for our z ~ 1 GRBs, leaving just a por- 
tion of the BB tail detectable in the XRT range. This may 
explain the lack of BB detection in GRBs 070419A and 
080319B, compounded by high absorbing column density 
in 070419A. Intrinsic column is likely to play an important 
role in BB detectability as discussed above and in Paper 
II. Power law photon index also shows a trend, with val- 
ues of r ~ 1.4 — 1.6 in the detected category while we find 
r ~ 2.1 — 2.9 in the possible BB category and something in 
between, F ~ 1.8 - 2.2, for the GRBs with no detected BB. 
Since the early power law index often evolves from hard to 
soft in GRBs (e.g. Zhang et al. 2007), it may be that the 
possible BBs are detected at later stages in the development 
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of the prompt emission. The steeper power law indices make 
determination of the slopes more difficult and typical errors 
on r are ± > 0.2 compared with ± ^ O.f for the detected BB 
category sources. The temporal power law index, however, 
has no bearing on the detection of a BB, with various light 
curve types sampled in each category (Figure [TJ. Consider- 
ing the prompt gamma-ray properties of these sources (as 
measured with Swift BAT for consistency across all GRBs, 
Sakamoto et al. 2011), while the first two BBs were discov- 
ered in subenergetic, long-duration GRBs, the candidates 
proposed since then have extended into the typical cosmo- 
logical GRB range. GRB 090618 is still the most energetic 
at _Eiso ~ 8 X 10®^ erg, and the span of _Eiso in the rest of 
our candidates is 10*^-10^^ erg. This strongly implies that 
thermal X-ray emission (or significant spectral curvature) 
can arise in any GRB, and its detection is dependent upon 
properties other than the prompt gamma-ray emission. 

We notice a distinctly lower BB luminosity and radius 
for the two low Ei^o GRBs 060218 and 100316D with de- 
tected BBs when compared with the GRBs with possible 
BBs, while GRB 101219B appears to lie in between the two 
categories. GRB 090618 does not follow the trend, with a 
very high BB luminosity, unsurprising given its very high 
afterglow fiux, but a relatively low inferred radius. The one 
BB parameter that correlates well with detection category is 
the fractional unabsorbed 0.3-10 keV X-ray flux contributed 
by the BB. This lies between 8 and 14% for the unambiguous 
cases, and between 32% and 65% in all possible detections. 
Interestingly, the restframe temperature of the BB remains 
similar for all detections, at around 0.2 keV, with the ex- 
ception of 090618 in which we see a somewhat higher kT, 
decreasing from 0.9-0.5 keV. 

For a number of GRB-SNe the SN properties have been 
measured (Hjorth & Bloom 2011; Cano et al. 2011b; Bufano 
et al. 2012; Sparre et al. 2011; Cobb et al. 2010; de Ugarte 
Postigo et al. 2011, in preparation; Richardson 2009; Ferrero 
et al. 2006; Melandri et al. 2012; Schulze et al. in prepara- 
tion), but beyond SN peak absolute magnitude the overlap 
with this sample is small since we also include photometric- 
only SN detections for which detailed ejecta information is 
not available. The GRBs with spectroscopically-determined 
supernovae are 060218, 081007, 091127, 100316D, 100418A, 
101219B and 120422A, spanning aU 3 categories of BB 
detection. SN peak absolute magnitudes are not available 
for 070419A, 081007 and 101219B. Host galaxy magnitude, 
which may be used as an approximate tracer for metallic- 
ity, is available for around half the sample: 060218, 060729, 
080319B, 090618, 091127 and 100316D hampering compar- 
ison with other properties. We can only say that 2 out 
of 3 GRBs with a clear BB have fainter hosts than the 3 
possible/no-BB GRBs, the bright host of 060218 being the 
exception. 

6.2 Origins of the thermal X-ray component 

The evolution of spectral parameters in our fits to 
GRB 101219B (Figure [S| may be compared with similar 
analyses performed for GRBs 060218 and 100316D (Figure 
4 of Starling et al. 2011; Figure 4 of Olivares et al. 2012; 
see also Kaneko et al. 2007) and GRB 090618 (Figure 8 of 
Page et al. 2011). While the passage of the spectral peak en- 
ergy through the X-ray band in the very soft GRBs 060218 



and 1001316D complicates the comparison, GRB 101219B 
shows a more similar behaviour to 090618 for most spectral 
parameters. The exception is BB luminosity, which steadily 
decreased in 090618 over 125-275 s since trigger. We find a 
general increase in 101219B from 160-544 s, followed by a 
decrease (though this decrease is only seen in the PC spec- 
tra which we stress can only be used as a guide to the BB 
behaviour). The increase does coincide with a fiattening of 
the light curve which may be interpreted as a short plateau, 
fiaring episode or other transition phase between the a ~ 1 
and a ^ 2 decays which may complicate the interpretation 
of our results. The luminosity variation has, however, little 
effect on the radius determination, which is almost exactly 
the same value and trend as is seen in 090618, increasing 
from OxlO'^^ cm to 7x10^^ cm in at most ~900s, (mean of 
556 s) implying an expansion velocity of 0.41q2C. Firstly, the 
radii that result are large compared with the 10^^ cm esti- 
mated for Wolf-Rayet stars (e.g. Cappa et al. 2004) , a likely 
progenitor type for long GRBs, but the existence of a thick 
wind could increase the radius probed in this way, as was 
suggested by some authors for the previous thermal X-ray 
components within the shock breakout model (but see e.g. 
Li 2007). Secondly, the inferred expansion velocity is very 
high, at least 20% of the speed of light, and together these 
considerations cast some doubt on the origin of this thermal 
emission in the shock breakout of GRB 101219B. We reiter- 
ate that these radii are illustrative only, since they assume 
both a single temperature BB and thermal equilibrium: at 
the high luminosities recovered here it seems likely a shock 
breakout would be highly relativistic, where thermal equilib- 
rium does not hold and indeed the predicted temperatures 
put the observable signature at energies above the soft X-ray 
band (Nakar & Sari 2012, before consideration of the effects 
of a wind). 

An alternative origin for thermal emission may be in 
a cocoon surrounding the jet. A number of our BB candi- 
dates are discovered during the steep decay phase of the 
X-ray light curve (Figure [T|. Pe'er, Meszaros & Rees (2006) 
showed that steep X-ray light curve decays (a > 2) can be 
produced by a relativistically expanding hot plasma cocoon 
associated with the jet. They demonstrate that some time 
after the prompt emission, typically a few hundred seconds, 
the bulk of the cocoon emission would fall in the X-ray band. 
The original spectrum is assumed to be a power law with 
power law index p = 2, as is expected, e.g., from synchrotron 
radiation. However, during the expansion phase of the co- 
coon the photons lose energy to the expanding plasma, and 
the spectrum is modified, resulting in a low energy ther- 
mal component combined with a power law over the 0.3- 
10 keV range. Origin in such a cocoon would allow measured 
BB radii to be larger than is expected for shock breakout, 
>10^^ cm, in line with many of our measurements. The high 
luminosities of the BB candidates would also be more easily 
explained in this model compared with the shock break- 
out model. The eventual predicted temporal X-ray decay is 
steep in this model, a ^ 2, which applies to some of our 
sample GRBs: 060729, 081007, 090618, 100418A, 101219B 
and 120422A. Interestingly, almost fiat early light curves are 
seen in GRBs 060218 and 100316D, i.e. in the most convinc- 
ing cases thus far for the shock breakout model. Addition- 
ally, very steep decays do not occur in the data in-hand for 
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070419A, 080319B or 091127, all GRBs where a BE could 
not be detected. 

Comparing the cocoon model of Pe'er et al. (2006) to 
our data for GRB 101219B we find that for reasonable model 
parameters both the overall spectral shape and the light 
curve decay can be reproduced. We have restricted the mod- 
elling to the time frame 100-250 s since trigger in this first 
approximate analysis, simply to see whether or not the gen- 
eral trends can be recovered. In Figure [5] we show some of 
these models, which stem from photons injected deep in the 
fiow having an initial photon energy distribution of a power 
law with index p — 2 above 10"^ me c^. The resulting spec- 
trum comprises low energy thermal spectra with modified 
power laws at higher energies, and the steep flux decay is 
produced, when the initial optical depth of the cocoon, r, is 
30-50. The Lorentz factor of the cocoon itself is found to be 
r = 4 — 7 with initial cocoon radius ~ 5 x lO'^^ cm~^. Total 
cocoon energy we have adopted, adjusted to match the light 
curve and energy of the early X-ray data, is of order a few 
X 10*^ erg (-Eiso as measured from prompt gamma-ray data 
is ~ 5 X lO''^ erg). 



7 CONCLUSIONS 

We have performed a systematic search for blackbody ther- 
mal components in the early Swift X-ray spectra of 11 GRBs 
with associated supernovae. This has resulted in the detec- 
tion of 4 compelling cases where, in addition to an absorbed 
power law, a BB-like component is required. We also note 
4 possible BB candidates and their caveats and 3 GRB-SN 
without detectable thermal emission. Amongst the detec- 
tions, we report for the first time details for GRB 101219B 
in which we see a cooling BB of restframe temperature 
~0.2 keV, luminosity ~3xl0''^ erg s~^ and inferred radius 
^4x10^^^ cm expanding at 0.2c or greater. Its general prop- 
erties appear to be intermediate between the low luminosity 
GRBs 060218 and 100316D and the classical GRB 090618, 
all of which have previously reported thermal X-ray emission 
that we also recover here. 

We cannot identify any strong correlations of BB prop- 
erties with the available prompt emission, SN or host prop- 
erties. We note some trends between BB 'category' and red- 
shift, power law spectral index and percentage contribution 
of the BB to the 0.3-10 keV flux; these likely do not refer to 
the intrinsic attributes of the GRBs, but rather to the de- 
tachability of the thermal component. All the possible BB 
candidates among our sample have large BB luminosities 
(> 10'*'^ erg cm"'^) and inferred radii (> 10^^ cm). These lu- 
minosities are difficult to reconcile with shock breakout mod- 
els since they imply highly relativistic outflows that would 
emit largely at energies above the Swift XRT bandpass. Our 
BB candidates occur during the steep decay phase of the X- 
ray light curve. We discuss a model for this emission arising 
in a relativistically expanding hot plasma cocoon associated 
with the jet. This model is able to generally reproduce the 
extreme BB parameters we retrieve from our fits, and is con- 
sistent with a steep temporal decay. Only GRBs 060218 and 
100316D do not show a temporal decay of a ^ 2 during the 
first ~1000s, where the BB luminosities and radii are signif- 
icantly lower than in the rest of our GRB-SN sample, and 
the shock breakout model has been argued by a number of 



Figure 6. Here we show models for emission from a relativisti- 
cally expanding hot plasma cocoon associated with the jet, de- 
scribed in Pe'er et al. (2006, see Section l6.2l l. These demonstrate 
that this cocoon model can produce the general shape seen in 
the spectrum (upper panel) and light curve (lower panel) of some 
of our GRB-SN thermal X-ray emission candidates. These are 
not fits to the data, but parameters are optimised for best rep- 
resentation of GRB101219B during Tq -I- 100 - 250 s. In each 
panel the three curves represent possible choices for model pa- 
rameters. Cyan: optical depth r = 30, initial expansion radius 
initial = 5 X lO^'^ Cm, total cocoon energy _Ec = 5 X 10'*® erg, 
r = 4.0; black: t = 30, nnitiai = 6 X lO^^ ^m, = 4.2 X 10*® 
erg, r = 4.5; red: t = 50, nnitial = 6 X lO^^ cm, = 4.2 X 10*® 
erg, r = 6.5. 




Observed energy [keV] 



10"' - 




10^ 10^ 

Observed time since burst [s] 



authors to be valid (Section [T}. Indeed, these two low lumi- 
nosity GRBs may be the product of quite different physics 
early on from the classical energetic GRBs that populate 
most of our sample, with GRB 101219B lying somewhere in 
the transition. 
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